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Abstract 

We define a percolation problem on the basis of spin configurations of the two dimensional XY 
model. Neighboring spins belong to the same percolation cluster if their orientations differ less than 
a certain threshold called the conducting angle. The percolation properties of this model are studied 
by means of Monte Carlo simulations and a finite-size scaling analysis. Our simulations show the 
existence of percolation transitions when the conducting angle is varied, and we determine the 
transition point for several values of the XY coupling. It appears that the critical behavior of this 
percolation model can be well described by the standard percolation theory. The critical exponents 
of the percolation transitions, as determined by finite-size scaling, agree with the universality class 
of the two-dimensional percolation model on a uniform substrate. This holds over the whole 
temperature range, even in the low-temperature phase where the XY substrate is critical in the 
sense that it displays algebraic decay of correlations. 

PACS numbers: 05.50.+q, 64.60.Cn, 64.60.Fr, 75.10.Hk 
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I. INTRODUCTION 



Consider the XY or planar model on the square lattice with periodic boundary conditions, 
described by the reduced Hamiltonian 

<ij> 

where the sum is over all nearest-neighbor pairs, and the Si are two-dimensional unit vectors 
labeled by the site number i. We restrict the nearest-neighbor interaction to be ferromag- 
netic, i.e., J > 0. 

When the temperature of this model is lowered, it undergoes a phase transition of an 



interesting character, which was explained by 
for the exponents were obtained by Nienhuis 



vosterlitz and Thouless [jj . More exact results 
for an 0(2) model in the same universality 



class. These results show that the renormalization exponent y t of the temperature is equal 
to 0, which means that the temperature-driven transition is of infinite order, i.e. the specific- 
heat singularity is extremely weak. In contrast, the magnetic susceptibility displays a very 
strong divergence when the temperature is lowered to the transition point. 

At temperatures T > below the transition point there is no spontaneous long-range 
order in the sense that the magnetization is zero 3fl. Instead, the low-temperature phase 
resembles a critical state; the correlations decay algebraically, with exponents that are still 
dependent on the temperature. 

Recently, there has been ample attention to percolation problems defined on a critical 
substrate, see e.g. Refs. 4j, 5, 6, 7 and references therein. It thus appears that, for Potts and 
O(n) models, the universal properties of such percolation transitions do reflect the nature of 
the critical substrate. These investigations are based on model representations with discrete 
degrees of freedom. 

It is thus interesting to investigate related problems using a substrate with continuous 
degrees of freedom. For instance, the mechanical properties of static granular matter can be 
analyzed in terms of the so-called force networks [a, |9|]. By introducing a threshold force, 
such that forces exceeding the threshold form clusters, a percolation transition is seen at a 
critical force threshold. This approach was applied to different models for granular piles, 
that are expected to belong to different universality classes. It was found that also the 



corresponding percolation behavior could discriminate between the different models 10] 



This might suggest that such substrate dependence is a general phenomenon for critical 
models with continuous degrees of freedom, i.e., the critical behavior of the percolation 
clusters might generally reflect the long-range correlations of the original degrees of freedom. 

In order to shed more light on this issue, the present work investigates a percolation prob- 
lem using the substrate of the two-dimensional XY-model. In particular we are interested 
how the universal properties of the percolation transition will depend on the temperature 
of the underlying XY model. We define the percolation problem such as to depend only 
on the XY configuration, and not on any additional random variables. This is achieved by 
introducing a "conducting angle" 9 such that neighboring spins whose orientations differ by 
less than 9 are connected by a percolation bond. This name is based on the analogy with 
a conductance problem of conducting units on a lattice, such that neighboring units are in 
electrical contact only if their orientations match to a sufficient degree. 

It is clear that this conducting angle also defines a threshold pair energy, below which a 
pair of neighboring XY spins is connected by a percolation bond. These percolation bonds 
define a bond percolation configuration involving a complete decomposition of the lattice in 
percolation clusters. An example of an XY configuration with the corresponding percolation 
cluster decomposition is shown in Fig. [TJ It is obvious that the resulting bond percolation 
configurations will not percolate for 9 = 0, and that they will percolate for 9 = ir and larger. 
We are thus left with the task to find the percolation threshold 9 C and the critical exponents 
as a function of the temperature. To answer these questions, we perform simulations of the 
two-dimensional (2D) XY spin model, which are described in Section [Til Section [TTT1 presents 
the numerical results, including the critical parameters. We include a short discussion in 
Sec. ED 



II. ALGORITHM 



For reasons of efficiency, we make use of a cluster algorithm lj| to simulate the XY 



12] to the model on 



model. We applied the single-cluster algorithm formulated by Wolff 
the square lattice. We recall the steps involved in one Wolff cluster flip: 

1. Choose an arbitrary direction as the y-direction, and denote by a, the angle between 
the spin si and the y-axis. Define an Ising spin Sj = ±1 with the same sign as the 
^/-component of s~l. Thus, the nearest-neighbor interaction term in Ji^ between spins 
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FIG. 1: Construction of a bond-percolation configuration (b) from an XY spin configuration (a). 
If the angle 9m between a pair (i, j) of neighboring spins is less than a given angle 9, the neighboring 
spins are connected by a bond, so that the lattice decomposes into a system of percolation clusters. 
The threshold was chosen as 9 = ir/6 in this figure. 

i and j reads 

Hij = -Ksi ■ Sj = —Ks*Sj - Ksjs^ = -Ks*Sj - K cos(«i) cos(<x,) Si-Sj , (2) 

where K = J/k&T. As far as the dependence of Tiij on the Ising variables Si and 
Sj is concerned, this is an Ising coupling between the two spins with strength Kij = 
K cos cVj cos a.,- . One may update the Ising variables using a cluster algorithm that 
takes into account these position-dependent couplings. The following steps are used 
to update the ^/-components of the spins. 

2. Choose a spin randomly, say on site i. For each nearest-neighbor site j of i, connect 
i and j by a bond with probability = max(0, 1 — e~ 2Kij ). Then do the same for 
each of the nearest-neighbor sites of each newly connected site, and so on. The process 
continues until no more new sites are connected. Then, the construction of the cluster, 
which contains all sites connected via some path of bonds to site i, is finished. 



3. Change the sign of the y-components of all spins in that cluster. 



III. NUMERICAL RESULTS AND ANALYSIS 



During the simulations, we constructed percolation cluster decompositions on the basis 
of a chosen conducting angle 9. For each such decomposition we sampled several quantities 
in order to estimate the second moment of the cluster size distribution S2, the probability P 
that there exists a cluster that wraps the system, a dimensionless ratio Q that can be related 
to the Binder cumulant, and the cluster size distribution function n s . We shall proceed to 
introduce these quantities in more detail and list their expected finite-size scaling properties. 

Since we are using finite systems with periodic boundary conditions, we define a "wrap- 



ping cluster" [13j as a percolation cluster that connects to itself along at least one of the 
periodic directions. For each XY configuration S we thus define a quantity p(S, 9) that has 
the value 1 if there exists a wrapping cluster, and p(S, 9) = otherwise. Thus, for a system 
with finite size L, the probability P(K, 9, L) that a wrapping cluster exists is given by 

P(K,9,L) = (p(S,9)}, (3) 

where the ensemble average is taken for an XY system with linear size L and coupling K. 
If there is a percolation transition at a conducting angle 9 C , one expects that 9 plays the 

n 

role of a temperature-like variable, and thus that the finite-size scaling behavior [14| of P is 
described by 

P(K, 9, L) = P (c) + a(9 - 9 c )L y v + ■■■ + b x L^ + b 2 L y * + • • • , (4) 

where y p controls the scaling of 9 and acts as a temperature-like exponent, corresponding 
with the bond dilution exponent in the language of the percolation model. The exponents y%, 
?/2, etc. are correction-to-scaling exponents, which are unknown in principle. The constant 
p( c \ which is defined as the value of P at 9 C in the limit L — > 00, and the exponent y p 
are universal, but the universality class of the present percolation problem remains to be 
determined. 

The second moment S2 of the percolation cluster size distribution can also be viewed as 
the mean size of the cluster containing an arbitrary point. It is defined as 

s ^w{t «?). < 5 > 

where Si is the size of the i-th cluster, N c the total number of clusters for a configuration, 
and N = L 2 is the volume of the system. The quantity 5*2 is also closely related with the 



generalization xrc of the (/-state Potts magnetic susceptibility to the random-cluster model, 
applied to the special case 5 = 1. This relation is expressed as xrc = NS2. 

In analogy with the Potts susceptibility, we expect the following finite-size-scaling behav- 
ior for S 2 in the neighborhood of a percolation threshold at 9 C : 

S 2 = L 2 ^- 2d [a + a x {9 - 9 c )Ifl> + ■■■ + hi* + b 2 I y * +•••], (6) 

where yh is the fractal dimension of critical percolation clusters. At the percolation threshold 
6> c , this equation reduces to 

S 2 = I 2 y h - 2d (a + 61^1 + b 2 V>* +••■)• (7) 



and Sa as 



We also define a dimensionless ratio related to the Binder cumulant 15] as 

Q = ^r — > ( 8 ) 

3^ 2) - 2S 4 

where is defined as 

The relation with the Binder cumulant is based on the fact that Eq. (jBJ) is obtained when 
the Binder ratio (m 2 ) 2 /(m 4 ) of magnetization moments of the Ising model is expressed in 
the language of the q = 2 random-cluster model. In analogy with the Binder cumulant, we 
expect that, in the neighborhood of the percolation threshold, it behaves like 

Q(K, 9, 1) = Qto + ai (0 - 9 C )U? + ■■■ + b x m + b 2 D>* +■■■ . (11) 



A. Wrapping probability as a function of 

We simulated XY systems with linear sizes 1 = 8, 16, 32, 64, 128, and 256 at different 
inverse temperatures: K = 0, 0.001, 0.01, 0.10, 5.00, and 10.00. For K = 0.9 and 
K = 1.0, system sizes up to L = 1024 and 2048 were simulated respectively. Typically, 
10 7 XY configurations were sampled. The samples were taken at intervals consisting of one 



0.9 1 1.1 1.2 1.3 1.4 1.5 

6 



FIG. 2: (color online). Wrapping probability versus 6 at K = 0.5. The curves are added as a guide 
for the eye. 

Metropolis sweep and four Wolff clusters. For each of these XY configurations, percolation 
cluster decompositions were constructed with several different values of the conducting angle. 

Some results for the wrapping probability P(K, 8, L) as a function of the conducting angle 
are shown in Figs. [2HH These figures show clear intersections corresponding with percolation 
thresholds. Furthermore, the behavior appears to be remarkably similar for different values 
of the XY coupling K. We applied the least-squares method to fit the wrapping probability 
by the finite-size-scaling equation fljj). The results for the percolation threshold 9 C are shown 
in Fig. [5] as a function of the XY coupling K. They are also included in Table [J for each 
value of K, together with the results for the exponent y p and the universal probability 
p( c \ The results for the latter quantity reproduce, within error bounds, the literature 
value P^ = 0.690473725 which applies to the ordinary two-dimensional percolation model 



16 



B. Dimensionless ratio Q as a function of 6 

We also sampled the ratio Q for systems at different inverse temperatures. Some of the 
data are shown in Figs. [6] to They display the same general behavior as the wrapping 
probability in the preceding subsection. 

A least-squares analysis of the data for Q(K,6,L) on the basis of the finite-size-scaling 
equation (TTTTl results in the estimates of y p , Q^ c \ and the percolation threshold 9 C . The 
estimates of y p and 6 C are in agreement with those found by fitting P(K, 6, L). The estimates 
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FIG. 3: (color online). Wrapping probability versus 9 at K = 1. The right-hand figure shows the 
details in the vicinity of 9 C , where the finite size effect becomes relatively strong. The curves are 
added as a guide for the eye. 
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FIG. 4: (color online). Wrapping probability versus 6 at K = 2.0. The curves are added as a guide 
for the eye. 



of the universal value are also listed in Table [T]for each K. Their values are in agreement 

0.87053 (2) for the ordinary two-dimensional percolation model. 



with the literature value 



C. The second moment S% at the percolation threshold 

We simulated the model with system sizes L — 8, 16, 32, 64, 128, 256 and 512 at the 
estimated percolation threshold 6 C for each K, and sampled the second moment S2 of the 
cluster-size distribution. Samples were taken at intervals consisting of one Metropolis sweep 
and four Wolff clusters. The numbers of samples taken for system sizes L = 8 to 32, L = 64 
to 128, and L = 256 to 512 are 10 7 , 4 x 10 6 and 2 x 10 6 respectively. For K = 1.0, which is 



10 

FIG. 5: Percolation threshold 9 C versus K. For K = 0, the critical value of the conducting 
angle is 1.61078 radians or about 92 degrees. This value applies to a system of randomly oriented 
conductors on the square lattice. The curve is added as a guide for the eye, and estimated error 
bars are smaller than the size of the symbols. 
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FIG. 6: (color online). Dimensionless ratio Q versus conducting angle 9 for K = 0.5. The curves 
are added as a guide for the eye. 

in the vicinity of the Kosterlitz-Thouless (KT) transition, additional simulations took place 
for system sizes L = 1024 and 2048, involving several times 10 6 samples per system size. 

Some of the data for S2 are shown as a function of the system size L in Figs. l9lfTT1 for 
some values of K. These figures use logarithmic scales, so that linear behavior means that 
5*2 behaves as a power of L, in accordance with criticality. For K = 1.0, which is close to 
the KT transitions, small deviations from linearity are visible. 

We fitted 62 by the finite-size-scaling formula Eq. ([7]), and thus obtained estimates of 
the fractal dimension yh, which are listed in Tabled for several values of K. For K = 1.0, 

n 
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FIG. 7: (color online). Dimensionless ratio Q versus conducting angle 9 for K = 1.0. The right- 
hand figure shows the details in the vicinity of 9 C . The curves are added as a guide for the eye. 
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FIG. 8: (color online). Dimensionless ratio Q versus conducting angle 9 for K = 2.0. The curves 
are added as a guide for the eye. 




FIG. 9: Second moment S% of the cluster size distribution versus finite size L, for XY coupling 
K = 0.5. The curve is added as a guide for the eye, and estimated error bars are smaller than the 
size of the symbols. 
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FIG. 10: Second moment £2 of the cluster size distribution versus finite size L, for XY coupling 
K = 1.0. The curve is added as a guide for the eye, and estimated error bars are smaller than the 
size of the symbols. 
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FIG. 11: Second moment S2 of the cluster size distribution versus finite size L, for XY coupling 
K = 2.0. The curve is added as a guide for the eye, and estimated error bars are smaller than the 
size of the symbols. 

a satisfactory fit could only be obtained by discarding system sizes L < 128 and including 
system sizes up to L = 2048. 

D. The distribution n s of the cluster size s 

At the percolation transition 9 C , the requirement that the cluster size distribution scales 
in a covariant way yields a power law for this distribution 

n s (s,9 c ) ocs- T , (12) 
1 1 
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FIG. 12: (color online). Cluster size distribution function as a function of cluster size s for K = 0.5, 
d c = 1.18179. 
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FIG. 13: (color online). Cluster size distribution function as a function of cluster size s for K = 1.0, 
6 C = 0.62210. 

where r is a critical exponent equal to 1 + d/yh, with d = 2 dimensions and the fractal 
dimension m of the percolation clusters. For the ordinary d = 2 percolation model, one has 
y h = 91/48 [19J, bo|, and thus r = 187/91. 

Parts of the data for n s are shown in Figs. n2lfT4l as a function of cluster size s, for 
different values of the XY coupling K. The axes in these figures use logarithmic scales. 
The data points are, in a wide range of cluster sizes, well approximated by straight lines, 
corresponding with a scale-invariant distribution according to a power law. The slopes 
of these lines are close to the value of r that applies to the ordinary percolation model. 
The curves in Figs. [T2lfT4l can obviously be brought to an approximate data collapse by 
introducing suitable prefactors, but it appears that there is still an appreciable finite-size 
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FIG. 14: (color online). Cluster size distribution function as a function of cluster size s for K = 2.0, 
9 C = 0.37560. 

dependence. 

For a numerical determination of r from the data for n s , we use a least-squares analysis of 
the range of s where n s is almost linear in Figs. rT2lfT4l thus excluding small clusters smaller 
than 20, and the largest clusters exceeding a size L Vh /10. Even then, two correction terms 
had to be included in order to obtain satisfactory residuals, according to the fit formula 



where the term with coefficient a 2 is significant only for K > 1. Satisfactory fits were ob- 
tained for all inverse temperatures. For example, for K = 0.75, we fitted the data simultane- 
ously for system sizes L = 32, 64,- • • , 512 using the above formula, and found r = 2.055 (2), 
yi = —0.57 (1) with x 2 P er degree of freedom almost equal to 1. For K = 2.0, a fit of the 
data for system sizes L = 32,64,- • • , 512 yielded r = 2.054 (2), y x = -1.0 (1), y 2 = -1.6 (1), 
again with a satisfactory value of x 2 as compared with the number of degrees of freedom. 

We thus obtained estimates of the exponent r that agree well with the value r = 187/91 = 
2.0549 • • • for the ordinary 2D percolation model, for several temperatures of the XY model. 

IV. DISCUSSION 

The numerical results presented in Sec. Ulllfor the exponents r, yt and y p , as well as for the 
universal probability and the universal ratio Q^ c \ agree accurately with the literature 
values applying to the two-dimensional percolation model. Thus our analysis shows the 




(13) 



TABLE I: Numerical results for the percolation threshold 9 C , the bond dilution exponent y p , the 
fractal dimension y^, the universal wrapping probability and the universal dimensionless ratio 
Q( c \ The exact values for two-dimensional percolation exponents are y p = 3/4 and yh = 91/48 = 
1.895833- ••. 
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existence of a transition in the ordinary percolation universality class, independent of the 
XY temperature. This result is as expected in the high-temperature phase where the XY 
spins display strong random disorder, but may seem somewhat surprising in the critical 
region including the low-temperature range. The problem discussed here can be viewed as 
a correlated percolation problem, in which the bond probabilities correlate as the nearest- 
neighbor differences in the XY model. For such long-range correlated percolation Weinrib 
2l| formulated a generalized Harris criterion to decide if the correlations are relevant for the 
percolation behavior or not. The critical behavior is expected to be in the universality class of 



ordinary (uncorrelated) percolation if au > 2 for correlations decaying with distance as r~ a , 
with v the percolation correlation length exponent v = 4/3. Indeed in the case at hand a = 2, 
and the correlations should thus be irrelevant. It stands in a strong contrast with a recent 
analysis [?| of a percolation problem defined on the basis of the two-dimensional O(n) model 
where n is a continuously variable parameter. For the 0(2) model, which belongs to the 
same universality class as the XY model, it was found [7J] that the percolation transition was 
driven by a bond dilution field that is only marginally relevant, i.e., y p = at that transition, 
while y p = 3/4 for ordinary percolation. It thus appears that the character of a percolation 
transition on a critical substrate depends on the precise definition of the percolation problem. 



The percolation problem of Ref. 7J was defined within regions separated by loops as defined 
in the context of the 0(n) loop model. These loops obviously display fractal properties js] 
at 0(n) criticality, and thus affect the nature of the percolation transition. In contrast, the 
spatial variation of the angle between neighboring XY spins is very smooth near criticality, 
and does not display boundary-like structures, even at the KT transition. While it is natural 
that the KT transition is reflected in some way in the behavior of the percolation transition 
at the critical value of the conducting angle, such effects are exposed by our numerical results 
only in terms of slow convergence in the analysis of the finite-size data. 
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